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We study the possibility of associating meta-stable Efimov trimers from three free Bose atoms in a tight trap 
realised, for instance, via an optical lattice site or a microchip. The suggested scheme for the production of these 
molecules is based on magnetically tunable Feshbach resonances and takes advantage of the Efimov effect in 
three-body energy spectra. Our predictions on the energy levels and wave functions of three pairwise interacting 
85 Rb atoms rely upon exact solutions of the Faddeev equations and include the tightly confining potential of an 
isotropic harmonic atom trap. The magnetic field dependence of these energy levels indicates that it is the lowest 
energetic Efimov trimer state that can be associated in an adiabatic sweep of the field strength. We show that the 
binding energies and spatial extents of the trimer molecules produced are comparable, in their magnitudes, to 
those of the associated diatomic Feshbach molecule. The three-body molecular state follows Efimov's scenario 
when the pairwise attraction of the atoms is strengthened by tuning the magnetic field strength. 
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I. INTRODUCTION 

Since the early days of quantum mechanics, the understand- 
ing of the complexity of few-body energy spectra has been the 
subject of numerous theoretical and experimental studies. Al- 
ready in 1935 Thomas Ql predicted that three particles may 
be rather tightly bound even when their short ranged pairwise 
interactions support only a single, arbitrarily weakly bound 
state. The energy of the tightly bound three-body state was 
found to diverge in the hypothetical limit of a zero range bi- 
nary potential. Thomas' discoveries were later generalised by 
Efimov 131, predicting that the number of bound states of three 
identical Bosons increases beyond all limits when the pairwise 
attraction between the particles is weakened in such a way that 
the only two-body bound state ceases to exist. Such an excited 
three-body energy state that appears under weakening of the 
attractive pairwise interaction is called an Efimov state. Con- 
versely, under strengthening of the attractive part of the binary 
potential an Efimov state disappears into the continuum. This 
remarkable quantum phenomenon of three-body energy spec- 
tra is usually referred to as Efimov's scenario. 

The existence of Efimov states in nature has still not been 
finally confirmed. Likely candidates may be found among the 
systems of identical Bosons, whose binary interaction poten- 
tial supports only a single, weakly bound state. Already in 
1977 Lim et al. |3] predicted the existence of an excited state 
of the helium trimer molecule 4 He3 which followed Efimov's 
scenario. This discovery has later been confirmed by inde- 
pendent theoretical studies 0,Hil0|l|l[l3[lll USi " g 
more accurate helium dimer potentials. From the experimen- 
tal viewpoint, the observations of Ref. fl^l clearly reveal that 
the helium dimer 4 He2 is indeed weakly bound. Even the he- 
lium trimer molecule has been detected by diffracting a helium 
molecular beam [13] from a micro-fabricated material trans- 
mission grating. While state selective diffraction experiments 
with helium trimers may, in principle, be possible l\M . there 
is still no conclusive evidence for the existence of their exited 



state. 

The possibility of manipulating the low energy inter-atomic 
interactions, using magnetically tunable Feshbach resonances, 
has provided new perspectives for the observation of Efimov's 
effect. Recent experiments with cold gases of Fermionic 
atoms as well as Bosonic species 

have demonstrated that adiabatic sweeps of the magnetic field 
strength can be used to associate highly excited diatomic 
Feshbach molecules with an arbitrarily weak bond. While all 
these experiments were performed in atom traps under com- 
paratively weak spatial confinement, there have been sugges- 
tions to produce molecules in the tightly confining light po- 
tential of an optical lattice l2lll . Since tight lattices suppress 
number fluctuations between different sites I22I1 . the molec- 
ular association may, in principle, be performed with two or 
even three atoms per site. A similarly tight or even stronger 
harmonic confinement of atoms may be achieved in microchip 
traps l23ll . The energy levels of a pair of interacting atoms in 
a tight micro-trap have been determined in Refs. 1241251 12611 . 
The universal properties l27t l28ll of three-body energy spec- 
tra in the presence of a confining harmonic potential have been 
studied in Ref. I29I1 . using an adiabatic approximation to solve 
the stationary Schodinger equation in hyper- spherical coordi- 
nates 03i • 

In this paper we exactly solve the Faddeev equations l3lll 
to determine the magnetic field dependence of the energy lev- 
els of three identical Bose atoms, whose pairwise interactions 
are tuned, using the technique of Feshbach resonances, in a 
harmonic micro-trap realised by an optical lattice site or a mi- 
crochip. Our results show that linear sweeps of the magnetic 
field strength can be used to populate the lowest energetic 
meta-stable Efimov trimer molecular state, largely in analogy 
to the association of diatomic Feshbach molecules. We show 
that the properties of the trimers produced, with respect to the 
strength of their bonds, are comparable to those of the associ- 
ated diatomic Feshbach molecule. We illustrate all our general 
findings for the example of the 155 G Feshbach resonance of 
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Rb Bose atoms. 

The paper is organised as follows: In Section[njwe discuss 
those universal properties of near resonant diatomic bound 
states that are crucial for the association of Efimov trimers. 
We then introduce the particular requirements on the type of 
Feshbach resonance, under which universality can be attained 
over a wide range of magnetic field strengths. Our results in- 
dicate that broad, entrance channel dominated Feshbach res- 
onances may be best suited to produce the trimers and pre- 
serve their stability on time scales sufficiently long to study 
their properties. We provide a general criterion that shows 
why the 155 G resonance of 85 Rb meets these requirements 

In Section [II]] we first illustrate the occurrence of the 
Thomas and Efimov effects in the energy spectra of three Bose 
atoms in free space, whose binary interactions are tuned using 
the technique of Feshbach resonances. Our discussion reveals, 
in particular, why all Efimov trimer states in such systems are, 
in general, intrinsically meta-stable. We then show how the 
three-body energy spectrum is modified due to the presence 
of the harmonic confining potential of an isotropic atom trap. 
These considerations allow us to identify the lowest energetic 
Efimov state as the molecular trimer state that can be associ- 
ated in an adiabatic sweep of the magnetic field strength. 

Section HVl illustrates the energy levels and wave functions 
of three 85 Rb atoms under the tight confinement of an optical 
lattice site or a microchip trap. We show that, under realis- 
tic conditions, the trimer molecules produced, when released 
from the lattice, are sufficiently confined in space that they 
can be identified as separate entities of a dilute gas. We then 
suggest a general scheme for their detection that directly takes 
advantage of the periodic nature of an optical lattice. 

All the details of our calculations are given in the appen- 
dices: Appendix 1X1 introduces the separable binary potential 
l35l l36l l37il that we have used to accurately describe the low 
energy scattering properties of a pair of 85 Rb atoms. We show, 
furthermore, how the separable potential approach can be ex- 
tended to determine the two-body energy levels, in the pres- 
ence of an isotropic harmonic atom trap, over a wide range 
of trap frequencies. We apply these techniques in Appendix 
151 to derive a general scheme to exactly solve the Faddeev 
equations for three pairwise interacting atoms including the 
confining trapping potential. Since our approach differs con- 
siderably from the known techniques for the solution of the 
Faddeev equations in free space, we provide a detailed de- 
scription of their numerical implementation. To demonstrate 
its predictive power with respect to three-body energy spectra, 
we provide estimates of the accuracy of our separable poten- 
tial approach through comparisons with ab initio calculations 
fioll of the helium trimer ground and excited state energies. 



II. ENERGY LEVELS OF A TRAPPED ATOM PAIR 

In this section we introduce the concept of universal two- 
body scattering and bound state properties characteristic for 
cold collision physics in the vicinity of zero energy reso- 
nances. These universal properties are crucial for the exis- 



tence of the Thomas and Efimov effects in three-body energy 
spectra. We describe the conditions under which the univer- 
sality of binary physical observables can be attained over a 
wide range of magnetic field strengths in experiments using 
Feshbach resonances to tune the inter-atomic interactions. We 
provide the relevant physical parameters of the microscopic 
binary potential that determine the energy spectra of an atom 
pair in free space as well as under the spatial confinement of 
an atom trap. We then show how the variation of the energies 
under adiabatic changes of the magnetic field strength can be 
used to associate diatomic molecules. Throughout this paper 
we discuss applications for the example of the 155 G Fesh- 
bach resonance of 85 Rb (1 G = 10~ 4 T). The underlying phys- 
ical concepts, however, are quite general and can be applied to 
cold collisions of other species of Bose atoms involving what 
we shall identify as entrance channel dominated resonances. 



A. Resonance enhanced scattering 

1. Magnetic field tunable Feshbach resonances 

Binary collisions in cold gases involve large de Broglie 
wavelengths which typically very much exceed all length 
scales set by the inter-atomic interactions. At the low col- 
lision energies the microscopic potential enters the descrip- 
tion of scattering phenomena only in terms of a single length 
scale, the s wave scattering length a. The experimental tech- 
nique of Feshbach resonances employs a homogeneous mag- 
netic field of strength B to manipulate the scattering length, 
taking advantage of the fact that the pairwise interaction de- 
pends on the coupling between the atomic Zeeman levels. 
Each Zeeman state is determined by the pair of total angular 
momentum quantum numbers (f, m /) of the hyperfine level 
with which the Zeeman state correlates adiabatically at zero 
magnetic field. In our applications to cold gases of 85 Rb the 
atoms are prepared in the magnetically trapped hyperfine state 
with the total angular momentum quantum number / = 2 and 
the orientation quantum number mj = —2 with respect to the 
direction of the magnetic field. The interaction between the 
atoms is usually described in terms of the binary scattering 
channels associated with the pairs of Zeeman levels of the in- 
dividual atoms. The relative energies between the dissociation 
thresholds associated with the channels can be tuned using the 
Zeeman effect. We shall denote the open s wave channel of 
a pair of asymptotically separated atoms of the gas as the en- 
trance channel. 

The typically weak inter-channel coupling can be grossly 
enhanced by tuning the energy E le!i (B) of a closed channel 
vibrational state |^ res ) (the Feshbach resonance level) in the 
vicinity of the dissociation threshold of the entrance chan- 
nel. General considerations l38l 13^1 show that a virtual en- 
ergy match between E les (B) and the threshold leads to a zero 
energy resonance in the entrance channel, i.e. a singularity of 
the scattering length, described by the formula: 



3 



Here «bg is usually referred to as the background scattering 
length, AB is the resonance width and B () is the position of the 
zero energy resonance. 



2. Universal properties of near resonant bound state wave 
functions 

The emergence of the zero energy resonance indicates the 
degeneracy of the binding energy £b(B) of the highest ex- 
cited vibrational multi-channel molecular bound state \if/^) 
(the Feshbach molecule) with the threshold for dissociation 
of the entrance channel at the magnetic field strength Bo- The 
magnetic field dependence of £b(B) is due to the perturba- 
tion of the physically relevant multi-channel energy levels 
by the strong coupling between the entrance channel and the 
closed channel Feshbach resonance state. The energy £t,(B) 
approaches the entrance channel dissociation threshold from 
the side of positive scattering lengths, while beyond the reso- 
nant field strength By, at negative scattering lengths, the bound 
state |^b) is transferred into a virtual state. At magnetic field 
strengths in the close vicinity of the zero energy resonance, 
all low energy binary collision properties are thus dominated 
by the properties of the bound state \t/fb), whose wave func- 
tion becomes universal in the limit a — > oo. This implies that 
the admixture of the Feshbach resonance level to the molecu- 
lar bound state \tf/\,) vanishes in accordance with the following 
asymptotic formula l39ll : 
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Here fi KS = dE res /dB is the virtually constant magnetic mo- 
ment of the resonance level |^ res ) and m is the atomic mass. 
We note that, in general, the physically relevant Feshbach 
molecular state \iff\,) and the Feshbach resonance level |0- res ) 
are considerably different with respect to their magnetic mo- 
ments and their spatial extents. The state \if/\,) binds the atoms, 
while l^res) may have only a short lifetime, in particular, when 
the inter-channel coupling is strong. 

In the vicinity of the magnetic field strength Bo the long 
range molecular bound state \i[r\,) consists mainly of its com- 
ponent in the entrance channel, which is given by the usual 
form of a near resonant bound state wave function 1 3 9L 14011 : 
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Its mean inter-atomic distance, i.e. the bond length, then di- 
verges like the scattering length in accordance with the rela- 
tionship: 



(r) 



J" d 3 rr\Mrt*a/2. 



(4) 



The associated binding energy is also determined solely in 
terms of the scattering length by the universal formula: 



Eb « -h 2 1 (ma 2 ). 



3. Entrance and closed channel dominated Feshbach resonances 

The possibility of magnetically tuning the scattering length 
using Feshbach resonances is unique among all physical sys- 
tems. The universal properties of the near resonant bound 
state wave function described by Eqs. (0, <E} and (0, how- 
ever, are not directly related to the inter-channel coupling and 
apply equally well, for instance, also to the deuteron in nu- 
clear physics lEUl and to the weakly bound helium dimer 4 He2 
van der Waals molecule Jl2ll . Since our applications to the 
association of Efimov trimer molecules crucially depend on 
the single channel nature of the bound state \tf>b) of the Fesh- 
bach molecule, we shall briefly outline the requirements on 
the Feshbach resonance that assure the validity of the univer- 
sal considerations over a significant range of magnetic field 
strengths. To this end, we consider the general properties 
of long range alkali van der Waals molecules that are deter- 
mined, to an excellent approximation, by the scattering length 
in addition to the asymptotic form -C(,/r 6 of the binary po- 
tential at large inter-atomic distances r. In accordance with 
Ref. El, we shall describe the dependence of the molecular 
energy level on the van der Waals dispersion coefficient C(, 
in terms of a mean scattering length a, which is given by the 
formula: 
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Here / v dw = 5 (inCe/h^j is usually referred to as the van 
der Waals length and F denotes Euler's gamma function. The 
binding energy of an alkali van der Waals molecule is then 
determined by the formula 14211 : 



Eb = -ft 2 1 |m(a — a) 2 J . 



(7) 



As discussed in detail for the examples of 23 Na and 85 Rb 
in Ref. 01, a variety of Feshbach resonances can be classi- 
fied on the basis of the properties of their associated Feshbach 
molecules: Throughout this paper, we shall denote a Feshbach 
resonance as entrance channel dominated when the binding 
energy Eb(B) is well approximated by Eq. Q in some range 
of magnetic field strengths about Bo, and, at the same time, 
the contribution of the mean scattering length a of Eq. im- 
proves the universal estimate of Eq. l|5}, i.e.: 



(8) 



The bound state \if/b) then describes a van der Waals molecule 
which implies that its properties are determined mainly by 
the entrance channel potential rather than the resonance level. 
In fact, the general considerations of Ref. ll39Tl show that the 
admixture of the closed channel resonance state |i^ res ) to the 
Feshbach molecular state \M) is negligible as soon as Eq. 
applies to its energy (cf., also, Fig. 3 of Ref. ll33lo . A gen- 
eral criterion for the applicability of Eq. Q can be derived 
from the two channel approach of Ref. 13911 which leads to the 
following inequality Il43ll : 
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Under the conditions of Eq. l|9}, the reason for the sup- 
pression of the closed channel contribution to the Feshbach 
molecule is the large detuning of the Feshbach resonance 
level from the dissociation threshold of the entrance chan- 
nel at the position Bo of the zero energy resonance 13^1 . 
The parameters for the 155 G Feshbach resonance of 85 Rb, 
i.e. UnJh = -3.12MHz/G H, C 6 = 4703 a.u. §1 (1 a.u. = 
0.095734 x 10- 24 Jnm 6 ), AB = 10.71 G 01, and a bg = 
-443 a (flo = 0.052918 nm) O, give the quantity in Eq. (|9} 
to be 4 x 10" 2 . The 155 G Feshbach resonance of 85 Rb is 
therefore entrance channel dominated. Similar conclusions 
were reached in Refs. 1I32H34I1 . 

There is also a variety of closed channel dominated Fesh- 
bach resonances, like, e.g., those in 23 Na l33ll . whose Fesh- 
bach molecular states |i/<b) become universal [cf. Eq. Q] only 
in a small region of magnetic field strengths about the zero 
energy resonance at Bo- These states are not intermediately 
transferred into a van der Waals molecule away from the reso- 
nance. Closed channel dominated resonances do not satisfy 
the criterion of Eq. and the bound states |^) are thus 
significantly influenced by the resonance level l^res) imme- 
diately outside the region of universality. Our considerations 
on three-body energy spectra take into account just a single 
scattering channel. We shall therefore focus in the following 
just on the entrance channel dominated Feshbach resonances 
of alkali Bose atoms. 



B. Adiabatic association of diatomic molecules in a tight atom 

trap 

Recent experiments have demonstrated the possibility of 
producing translationally cold diatomic Feshbach molecules 
using linear downward ramps of the Feshbach resonance 
level across the dissociation threshold of the colliding atoms. 
These studies on molecular association have been performed 
in quantum degenerate Fermi gases fl5[ fl6ll . consisting of 
an incoherent mixture of two spin states, as well as in dilute 
vapours of cold Bose atoms 1 1 7[ 1 1 8l 1 1 9l 12011 . 

The highly excited Feshbach molecules produced can be 
quite unstable with respect to de-excitation upon collisions 
with surrounding atoms. Exact calculations of the de- 
excitation rate constants are challenging and have been per- 
formed only for transitions between tightly bound states 14711 . 
Most of the current knowledge, therefore, relies upon experi- 
mental evidence. For Fermionic species it has been predicted 
|j48ll that the de-excitation mechanism is particularly efficient 
when the bond length of the Feshbach molecule is sufficiently 
small for its wave function to have a significant spatial overlap 
with more tightly bound molecular states. The experimental 
studies of Ref. 14911 confirm this trend. 

The observation of large collisional de-excitation rate con- 
stants of Feshbach molecules consisting of Bose atoms has 
been reported in Ref. ll50ll . The 23 Na resonance studied in 
these experiments, however, is closed channel dominated 1 33]. 
The properties of the 23 Na2 Feshbach molecules of Ref. 15011 
are thus rather different from those that we shall discuss in the 
following applications l33ll . The general experimental trends 



for both the Fermionic Jl5ll and the Bosonic l20ll species sug- 
gest that broad, entrance channel dominated Feshbach reso- 
nances may be best suited to associate a large portion of the 
atoms to Feshbach molecules and stabilise them. Most of the 
currently known entrance channel dominated Feshbach reso- 
nances have been found, however, in Fermionic gases Jl5ll5lll . 

Inelastic de-excitation collisions with background atoms 
may be efficiently suppressed when the molecules are pro- 
duced in the tight micro-traps of an optical lattice with an av- 
erage occupation of two or three atoms per site, respectively. 
Since the production of diatomic Feshbach molecules and Efi- 
mov trimers in tight harmonic atom traps can be performed in 
similar manners, we shall first discuss the underlying physi- 
cal concept for the simpler case of the association of a pair of 
atoms. 

The association of diatomic Feshbach molecules using 
magnetically tunable Feshbach resonances is closely related 
to the variation of the two-body energy spectrum under adia- 
batic changes of the magnetic field strength. Figure [2 shows 
such an energy spectrum versus the magnetic field strength 
B for a pair of 85 Rb atoms in a tight spherically symmetric 
harmonic atom trap, whose high frequency of v bo = 300 kHz 
may be realised in a microchip trap 15211 or in an optical lattice 
12511 . The details of the underlying calculations are explained 
in Appendix 1X1 In Fig.[2the binding energy of the Feshbach 
molecule in free space is shown for comparison. The spherical 
symmetry of the atom trap allows us to separate the centre of 
mass from the relative coordinates of the atoms ll24ll25ll26ll . 
and only the relevant levels of the relative motion are depicted 
in Fig. [2 We have chosen the zero of energy, for each mag- 
netic field strength, at the dissociation threshold of the en- 
trance channel in free space. Following the adiabatic curves 
of the energies clearly reveals that a pair of 85 Rb atoms occu- 
pying the level closest to the threshold on the low field side of 
the zero energy resonance is transferred into the level of the 
Feshbach molecule when the magnetic field strength is varied 
adiabatically across the resonance. The energy of the trapped 
molecular level and the free space binding energy E\, approach 
one another as the magnetic field strength is increased. A sim- 
ilar statement applies to their wave functions. This reflects 
the physical concept of the adiabatic association of diatomic 
molecules in tight atom traps. 

The two-body trap level closest to the dissociation threshold 
can, in principle, be prepared by loading a Bose-Einstein con- 
densate adiabatically into an optical lattice. In the case of 85 Rb 
the negative background scattering length of Ob g = -443 «o 
prevents such a condensate from being stable on the low field 
side of the resonance. To produce 85 Rb2 Feshbach molecules 
via an adiabatic sweep of the magnetic field strength, the 
Bose-Einstein condensate needs to be prepared on the high 
field side of the resonance before it is loaded into the lattice. 
The resonance then needs to be crossed as quickly as possible 
to reach its low field side and, at the same time, avoid a sig- 
nificant heating of the atomic cloud l33ll . Such a sequence of 
magnetic field sweeps is described in Ref. l20ll in the context 
of the adiabatic association of 85 Rb 2 Feshbach molecules in a 
cold gas. 
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FIG. 1 : (Color online) Magnetic field dependence of the vibrational 
energy levels E 2 b of a pair of 85 Rb atoms in a v ho = 300 kHz trap (up- 
per part) as compared to the binding energies of the highest excited 
vibrational state of the 85 Rbi dimer molecule in free space (lower 
part). The circles in the upper part indicate numerical solutions of the 
two-body Schrodinger equation for a microscopic interaction poten- 
tial explicitly incorporating the exact scattering length as well as the 
exact asymptotic -C^/r 6 interaction energy, while the solid curves 
indicate calculations using the separable potential approach of Ap- 
pendix^ The diamonds in the lower part are experimental binding 
energies obtained from Ref. 14611 and the squares indicate results of 
full coupled channels calculations of S. Kokkelmans l44ll . The solid 
curve indicates the dimer binding energies obtained from the separa- 
ble potential approach, while the dashed dotted curve corresponds to 
their near resonant approximation of Eq. The inset of the lower 
part of the figure shows the singularity of the scattering length at the 
magnetic field strength B = 155.041 G I4H1 . 



III. THOMAS AND EFIMOV EFFECTS 

In this section we discuss the Thomas and Efimov effects in 
the three-body energy spectra of interacting Bose atoms. We 
show that these phenomena occur when the binary scattering 
length is tuned by a magnetic field in the vicinity of a zero 
energy resonance. We then describe how the Efimov spectrum 
is modified in the presence of a trapping potential. Our results 
indicate that it is the lowest energetic Efimov trimer state that 
can be populated by adiabatic changes of the magnetic field 
strength. 



A. Three-body energy levels in free space 

Throughout this section, we consider three Bose atoms that 
interact pairwise through their binary potential. This assump- 
tion is justified for the description of weakly bound molecules, 
like the Efimov trimers in the present applications, whose 
inter-atomic separations very much exceed the van der Waals 
length. Tightly bound trimer molecules may be significantly 
influenced by genuinely three-body forces, which we shall ne- 
glect in the following. We assume furthermore that the binary 
potential supports at most a single bound state, the Feshbach 
molecule \iff\,). We thus neglect all tightly bound diatomic 
states, whose spatial extents are typically much smaller than 
the van der Waals length. In view of the large separation of 
the length scales between Efimov trimers and the Feshbach 
molecule on the one hand and the tightly bound dimer states 
on the other hand, we believe that this treatment provides an 
excellent approximation to the three-body states and their low 
energies we consider in this paper. We note, however, that any 
trimer state with an energy above the two-body ground level 
can, at least in principle, decay into a dimer bound state and a 
third free atom in accordance with energy conservation. As, in 
free space, the binding energies of three atoms are thus strictly 
limited, from above, by the two-body ground state energy, the 
trimer molecules under consideration are all in meta-stable 
states. Their associated lifetimes may depend sensitively on 
the details of the binary and three-body interactions. 

Given that the binary interactions support just a single, ar- 
bitrarily weakly bound state, it has been predicted, in terms 
of a rigorous variational treatment by Thomas I 1 !], that three 
particles can be comparatively tightly bound. The three-body 
ground state can persist even when the binary interactions are 
weakened in such a way that their only bound state ceases to 
exist. Such three-body bound states that exist in the absence 
of any bound two-body subsystem are usually referred to as 
Borromean states 15 311 . 

The Thomas scenario of Borromean states has been sub- 
sequently generalised by Efimov 01, in a striking way, pre- 
dicting that the number of three-body bound states of identi- 
cal Bosons increases beyond all limits when the energy of the 
only two-body bound state is tuned towards the dissociation 
threshold. Efimov's effect is closely related to the spatial ex- 
tent of the near resonant two-body bound state wave function, 
which is determined by the scattering length in accordance 
with Eq. 0. According to Efimov's treatment, it is indeed 
the scattering length, rather than the range of the potential, 
that sets the scale of the range of the effective three-particle 
interactions at the low collision energies under consideration. 
When the two-body binding energy reaches the dissociation 
threshold these interactions therefore acquire a long range. In 
contrast to a short range binary potential, however, long range 
interactions can support infinitely many bound states 15411 . 
All these Efimov states are spherically symmetric and their 
energies accumulate at the three-body dissociation threshold. 
Their number is predicted to follow, in the limit \a\ — > 00, the 
asymptotic relationship: 

AWov £ - log (Pc \a\ IK) . (10) 

71 
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Here p c is a momentum parameter related to the range of the 
binary interactions. Efimov's remarkable results have been 
subsequently confirmed by Amado and Noble l55ll . 
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FIG. 2: (Color online) Magnetic field dependence of the vibrational 
energy levels of 85 Rb3 trimers relative to the binding energy Et, of 
the 85 Rb2 Feshbach molecule (using E\, = for B < Bo) in free 
space (cf. Fig.0. The first Efimov states, whose energies are indi- 
cated by the solid and dashed curves, emerge at about 154.4 G and 
155 G, respectively. The energies of the other Efimov states are not 
resolved even on the logarithmic energy scale. The second Efimov 
state (dashed curve) ceases to exist at about 155.5 G. 



dissociation threshold, i.e. the binding energy E\, of the Fesh- 
bach molecule. The energy £3 of the second Efimov state thus 
crosses the dimer binding energy at a magnetic field strength 
of 155.5 G. 

Following the rigorous proof of Efimov's effect by Amado 
and Noble 15511 . the parameter p c of Eq. i\0\ may be esti- 
mated, using the separable potential approach to the low en- 
ergy spectrum of a pair of alkali atoms of Appendix|X] to be: 



p c « 2%/(na). 



(11) 



In agreement with Thomas' 1 1 ] and Efimov's [2] original sug- 
gestions, Eq. (II It recovers the order of magnitude of /?/r e ff, 
where r eff = ±[r(l/4)/r(3/4)] 2 a * 2.9 x a is the effective 
range of the interaction between a pair of alkali atoms in the 
limit of large scattering lengths 1571 15811 . Equation dl 11 1 thus 
confirms that, in contrast to the associated two-body prob- 
lem, the low energy physics of three Bose atoms crucially 
depends on the range of the binary potential. In fact, in the 
hypothetical limit of a zero range potential not only the num- 
ber of three-body bound states becomes infinite but also the 
three-body ground state energy diverges y|. This singular be- 
haviour clearly reveals that the low energy three Boson prob- 
lem is unsuited for a treatment in terms of pairwise contact 
interactions in the absence of energy cutoffs. 



B. Adiabatic association of Efimov trimers in an atom trap 



Efimov's scenario can be realised, using the technique of 
Feshbach resonances, by magnetically tuning the binary scat- 
tering length of three Bose atoms in the vicinity of a zero en- 
ergy resonance. Figure [2] shows the energy levels of three 
85 Rb atoms in free space versus the magnetic field strength 
in the vicinity of the 155 G Feshbach resonance. The exact 
three-body binding energies have been determined using the 
momentum space Faddeev approach l56ll and the separable 
binary potential of Appendix [X] According to Eq. fllOi , the 
number of Borromean Efimov states at negative binary scat- 
tering lengths increases beyond all limits when the pairwise 
attraction is strengthened in such a way that the two-body 
bound state emerges at the dissociation threshold. Two of 
their energies, i.e. the solid and dashed curves, are resolved 
on the logarithmic scale of Fig. [2 The dotted dashed curve is 
associated with the energy of the comparatively tightly bound 
Borromean state predicted by Thomas QJ on the low field side 
of the zero energy resonance at Bo = 155.041 G (vertical solid 
line). 

It turns out that, as the bond of the dimer state is strength- 
ened any further, the energies of the Efimov states succes- 
sively cross the two-body binding energy and become un- 
bound. Beyond the crossing point, the Efimov states can 
decay into a bound two-body subsystem, i.e. the Feshbach 
molecule, and a free particle, in accordance with energy con- 
servation. This explains why the number of three-body bound 
states decreases as the attractive pairwise interactions are 
strengthened in the presence of a two-body bound state. In 
Fig. [2] we have chosen the zero of energy to be the three-body 



The spatial confinement of an atom trap restricts the bond 
length of the Feshbach molecule. This, in turn, implies that 
the energy levels of three trapped atoms, unlike those in free 
space, do not have any accumulation point even when the 
magnetic field strength is tuned across a singularity of the 
binary scattering length. Figure reveals, however, that the 
energy spectrum of three 85 Rb atoms depends sensitively on 
the magnetic field strength, largely in analogy to the two-body 
spectrum of Fig.^ The energy levels in Fig.[3]have been ob- 
tained from exact solutions of the Faddeev equations, in the 
presence of a spherically symmetric trapping potential, using 
the approach of AppendixlBl 

We have chosen the low frequency Vho = 200 Hz of a typ- 
ical magnetic atom trap in Fig.|3j which allows us to directly 
compare the energy levels of three 85 Rb atoms in the pres- 
ence of the trap with those of the Efimov states in Fig. [2 This 
comparison shows that the trapped three-body energy level En 
(solid curve), closest to the three-body dissociation threshold 
in free space (dotted horizontal line), correlates adiabatically, 
in the limit Vho — » 0, with the energy of the first Efimov state. 
A magnetic field pulse sequence similar to the one discussed 
in Subsection lH Bl and in Ref. 1 20] can, in principle, be used to 
populate the trapped three-body energy level associated with 
the solid curve in Fig.[3]on the low field side of Bo. An adia- 
batic upward sweep of the magnetic field strength across the 
three-body zero energy resonance of Fig. [2] at about 154.4 G 
then transfers the trapped state into the first Efimov state. The 
physical concept underlying the adiabatic association of Efi- 
mov trimers in atom traps is, therefore, completely analogous 
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FIG. 3: (Color online) Magnetic field dependence of the vibrational 
energy levels £33 of three interacting 85 Rb atoms in a low frequency 
Vh = 200 Hz trap as compared to the energies of the free space Efi- 
mov states of Fig. [5] The energies of the trapped atoms are shown 
relative to E2B + f^Vho. where we have chosen E2B as the lowest 
energy of a pair of trapped 85 Rb atoms, i.e. the two-body level that 
correlates adiabatically, in the limit Vh — > 0, with the binding energy 
E h of the Feshbach molecule of Fig. Q In analogy, we have sub- 
tracted the dimer binding energy E b from the Efimov trimer energies 
in free space for the purpose of comparison. 



trends in their dependence on the magnetic field strength. The 
variation of the trap frequencies mainly affects the spacings 
between the energy levels. In fact, all our considerations with 
respect to the association of Efimov trimer molecules depend 
just on the possibility of trapping exactly three atoms rather 
than on the tightness of the spatial confinement. 
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FIG. 4: (Color online) Magnetic field dependence of the vibrational 
energy levels £ 3B of three 85 Rb atoms relative to the zero point en- 
ergy jhvha of three hypothetically non-interacting atoms in a tightly 
confining Vho = 300 kHz trap. 



to the considerations on the production of diatomic Feshbach 
molecules in Subsection lllBI 

IV. TRIMER MOLECULES IN AN OPTICAL LATTICE 

In this section we describe three interacting 85 Rb atoms in 
a tight micro-trap of an optical lattice site or in a microchip 
trap with a realistic oscillator frequency. Our results indicate 
that the spatial extent of the wave function of the first Efimov 
state can be tuned in such a way that it is smaller than realistic 
mean inter-atomic separations of dilute gases. We discuss, 
furthermore, how the periodicity of an optical lattice can, in 
principle, be used to detect the trimer molecules. 

A. Three-body energy levels and wave functions in tightly 
confining atom traps 

The sites of an optical lattice, in general, confine the atoms 
much more tightly than usual magnetic traps for cold gases. 
Their harmonic frequencies differ by several orders of magni- 
tude, from tens to hundreds of kHz in the case of a lattice 1251 
or a microchip trap 12 3l 15211 as compared to about 100 Hz for a 
conventional magnetic trap. During the course of our studies, 
we have calculated the energy levels of three 85 Rb atoms for a 
variety of trap frequencies extending from 200 Hz in Fig.|3]to 
1 MHz. Our results indicate that, despite the pronounced dif- 
ferences in the spatial confinement, all spectra follow the same 



As a typical example of our results, Fig. @] shows the en- 
ergy levels of three 85 Rb atoms versus the magnetic field 
strength in a tightly confining atom trap with a frequency of 
= 300kHz. We note that in Fig.0]the energies are given 
relative to the zero point energy of three hypothetically non- 
interacting trapped atoms, while in Figs.[2]and[3]we have cho- 
sen the zero of energy at the three-body dissociation threshold 
in free space, i.e. at the binding energy E\, of the Feshbach 
molecule. The solid curve of E2 in Fig.|4]thus correlates adi- 
abatically, in the limit Vho — » 0, with the first Efimov state of 
Figs. [2]and[3] Figure reveals that the trapped first Efimov 
state is transferred into a meta-stable trimer molecule, within 
a range of magnetic field strengths from 154.4 G to at least 
160 G, when it is adiabatically released from the trap. The 
modulus of its binding energy is always slightly larger than 
|i£bl- These predictions suggest that the adiabatic association 
of the first Efimov trimer state is, in principle, feasible and 
leads to reasonably strong bonds. We note that the produc- 
tion of the weakly bound 85 Rb 2 Feshbach molecules has been 
observed over a wide range of magnetic field strengths ll2oTl . 

Weakly bound molecules, such as Efimov trimers, are char- 
acterised by a large spatial extent of their wave functions iB . 
In order to identify them as separate entities in a dilute gas, it 
is crucial for the size of the molecules to be much smaller than 
the mean spacing between their centres of mass. As discussed 
in Appendix |EJ even the isotropic wave functions of the Efi- 
mov trimers depend on three parameters and can therefore not 
be directly visualised. To give an impression of the spatial 
structure of their wave functions, Fig.|5]shows the hyper-radial 
probability densities P(R) of the three lowest energetic 85 Rb 
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FIG. 5: (Color online) Hyper-radial probability densities P(R) of the 
lowest energetic vibrational states of three 85 Rb atoms at 6 = 158. 1 G 
for the trap frequencies Vh = 1 MHz, v^ = 300 kHz, and Vh = 
50 kHz. The legends show the associated energies of the interacting 
atoms relative to the zero point energy of three hypothetically non- 
interacting atoms, i.e. E 3B - ^hvho- The hyper-radius is given on a 
logarithmic scale. 



trimer states, at a magnetic field strength of 158.1 G, for trap 
frequencies decreasing from 1 MHz down to 50 kHz. In anal- 
ogy to the case of diatomic molecules, the number of zeros 
of P(R) is related to the degree of excitation of the three-body 
energy state. The solid curves are associated with the first 
Efimov state, which can be produced through an adiabatic up- 
ward sweep of the magnetic field strength. The wave func- 
tion of this state extends over a few hundreds of Bohr radii in 
the 1 MHz trap and completely decays at about 2000 a.u. in 
the case of the 50 kHz trap. The sizes of the trimer states 
produced are thus comparable to the single particle oscillator 
lengths (a ho = [h/(mv ho )] l/2 a 1000 a.u. for 85 Rb atoms in 
a 50 kHz trap). This indicates that the three atoms occupy a 
volume comparable to the size of the single atom oscillator 
ground state. Although these length scales very much exceed 
the spatial extent of even the most loosely bound observed 
diatomic ground state molecule (the helium dimer fcHo . re- 



alistic cold gases are usually sufficiently dilute that the wave 
functions of these Efimov trimers do not overlap each other 
once they are adiabatically released from the lattice (cf., e.g., 
Ref. 15^1 for an estimate of the remarkably large bond lengths 
of 6 Li2 Feshbach molecules produced in a dilute gas). 



B. Detection of Efimov trimer molecules in an optical lattice 

A variety of present day detection techniques for weakly 
bound Feshbach molecules in cold gases relies upon direct rf 
photo-dissociation spectroscopy 1151 . atom loss and recovery 
measurements 

El 

12011 . or the spatial separation of the molec- 
ular cloud from the remnant atomic gas 1171 1 1 8t 1 1 911 . While 
all these techniques may be applicable, in one way or another, 
also to Efimov trimers, an optical lattice lends itself to a rather 
different approach to their detection; the mass spectrometry 
using the periodic light potential as a diffracting device. We 
shall focus in the following on the perspectives of the diffrac- 
tion technique. 

The crucial coherence properties of Bose-Einstein con- 
densed atomic gases loaded into optical lattices have been 
studied in detail in Ref. l22ll . These experiments indicate that 
the possibility of diffraction depends on the weakness of the 
light potential. Contrary to this, the association of Efimov 
trimers requires high tunnelling barriers between the individ- 
ual sites to protect the molecules against inelastic collisions. 
In accordance with Ref. I22I1 . we expect that adiabatically re- 
leasing the lattice depth transfers the gas of the Efimov trimers 
produced from their insulating phase back into the super-fluid 
phase. Super-fluid gases in optical lattices, however, can be 
diffracted 12211 . 

The spatial periodicity of the optical lattice implies that the 
momenta of the cold molecules, when released from the light 
potential, are determined by multiples of the reciprocal lat- 
tice vectors with a negligible spread under the conditions of 
super-fluidity. The associated quantised velocity transfers ob- 
tained from the lattice, and consequently also the diffraction 
angles, are therefore inversely proportional to the molecular 
mass. The principle of this mass selection technique has been 
demonstrated in several earlier experiments on the spatial sep- 
aration of weakly bound helium dimers fl2l fl3ll and trimers 
fl3ll as well as sodium dimers Ir30ll from molecular beams. 
The freely expanding 85 Rb Efimov trimers of the present ap- 
plications may then be dissociated and imaged using the meth- 
ods of Refs. II6II17III8LI19LI20II . The mass selection detection 
technique suggested in this paper is general and could be ap- 
plied, for instance, also to diatomic Feshbach molecules. 

We expect that the most serious constraint on the produc- 
tion and detection of Efimov trimers in optical lattices is their 
intrinsic meta-stability. In the special case of 85 Rb 3 molecules 
there are two mechanisms that can lead to their spontaneous 
dissociation: The first mechanism consists in the decay into a 
fast, tightly bound dimer molecule and a fast atom. The sec- 
ond decay scenario involves spin relaxation of a constituent of 
the trimer molecule in analogy to the studies of Refs. l2(ilr)lll . 
While the present considerations do not allow us to estimate 
the molecular lifetimes associated with the first decay mecha- 
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nism, it has been shown in Refs. that spin relaxation 

can be efficiently suppressed by increasing the spatial extent 
of the molecules. It can be ruled out completely for other 
atomic species, in which the individual atoms are prepared in 
their electronic ground state. We are, however, not aware of 
other well studied entrance channel dominated Feshbach res- 
onances of identical Bose atoms besides 85 Rb. 



V. CONCLUSIONS 

We have studied in this paper the association of weakly 
bound meta-stable trimer molecules from three free Bose 
atoms in the ground state of a tight micro-trap of an optical lat- 
tice site or of a microchip. Our approach takes advantage of a 
remarkable quantum phenomenon in three-body energy spec- 
tra, known as Efimov's effect. Efimov's effect occurs when the 
binary scattering length is tuned in the vicinity of a zero en- 
ergy resonance and involves the emergence of infinitely many 
three-body molecular states. We have shown that this scenario 
can be realised by magnetically tuning the inter-atomic inter- 
actions using the technique of Feshbach resonances. 

The association scheme for trimer molecules, suggested in 
this paper, involves an adiabatic sweep of the magnetic field 
strength across a three-body zero energy resonance and can 
be performed largely in analogy to the well known associa- 
tion of diatomic molecules via magnetically tunable Feshbach 
resonances. Our results indicate that the predicted binding 
energies and spatial extents of the Efimov trimer molecules 
produced are comparable in their magnitudes to the associ- 
ated quantities of the diatomic Feshbach molecules. We have 
illustrated our general considerations for the example of 85 Rb 
including a suggestion for a complete experimental scenario 
and a possible detection scheme. 

Once the meta-stable trimer molecules are produced, the 
possibility of tuning the inter-atomic interactions using Fesh- 
bach resonances may, in principle, be exploited to study the 
Efimov property of the trimer state. Since, according to the 
predictions of this paper, it is the first Efimov state that gets as- 
sociated in an adiabatic sweep of the magnetic field strength, 
the trimer molecules are expected to dissociate as the pair- 
wise attraction between the atoms is strengthened. This, at 
first sight, counterintuitive scenario can be realised by tuning 
the magnetic field strength away from the zero energy res- 
onance on the side of positive scattering lengths. Once the 
energy of the diatomic Feshbach molecule crosses the bind- 
ing energy of the trimers, the Efimov states dissociate into 
the Feshbach molecule as well as a third free atom. In this 
way, the technique of Feshbach resonances could provide a 
unique opportunity to finally confirm this predicted |2] but as 
yet unobserved fascinating quantum phenomenon in the en- 
ergy spectrum of three Bose particles. 



APPENDIX A: SEPARABLE TWO-BODY INTERACTION 

In this appendix we provide a convenient effective binary 
interaction potential that well describes the relevant low en- 
ergy vibrational levels of an atom pair, both in free space and 
under the strong spatial confinement of a tight micro-trap of 
an optical lattice site or of a microchip. We determine the 
parameters of the effective potential in terms of the s wave 
binary scattering length a and the van der Waals dispersion 
coefficient C(,, which characterises the interaction energy at 
asymptotically large inter-atomic distances. 



1. Overview of the separable potential approach in free space 

a. Hamiltonian 

We first consider a pair of identical Bose atoms of mass m 
at the positions ri and r 2 , respectively, which interact via the 
microscopic potential V in the absence of a confining atom 
trap. The associated free-space binary Hamiltonian is then 
given by the formula: 



#free - #o' ee + ^ 

Here the non-interacting Hamiltonian 



wr = - 



n 2 



2(2m) 



■vi- 



2(f) 



(Al) 



(A2) 



accounts for the kinetic energy, where R = ^(rj + r{) and 
r = ri -T2 denote the centre of mass and relative coordinates 
of the atom pair, respectively. We assume in the following 
that the centre of mass is at rest and focus only on the relative 
motion of the atom pair. The non-interacting Hamiltonian in 
free space then reduces to H^' ee = —fi 2 V^/m. 



b. Transition matrix 

In our subsequent applications to three-body systems it will 
be convenient to represent all bound and free energy levels of 
the binary subsystems in terms of their transition matrix T(z) 
associated with the relative motion of the atoms, whose sin- 
gularities as a function of the continuous variable z determine 
the two-body energy spectrum. In general, the transition ma- 
trix associated with the interaction potential V can be obtained 
from the Lippmann-Schwinger equation (62]: 

(A3) 



f(z) = V + VG f ™(z)f (z). 



Here Go" ee (z) = (z - H^y 1 is the Green's function of the 
relative motion of the atoms in the absence of inter-atomic 
interactions. 
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The solution of the two-body Lippmann-Schwinger equa- 
tion ( IA3I for the microscopic inter-atomic potential V(r) is 
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a demanding problem in its own right. The full binary in- 
teraction, however, describes a range of energies much larger 
than those accessible to cold collision physics. We shall there- 
fore introduce a simpler effective potential which properly ac- 
counts just for the relevant low energy physical observables. 
Considerations 13^1 beyond the scope of this paper show that 
all physical observables associated with cold binary collisions 
can be described by a single parameter, the s wave scattering 
length a. The T matrix determines the scattering length by its 
plane wave matrix elements in the limit of zero energy: 



<p' = 0\f (0)|p = 0) = 



1 Anh 2 



{2nKy m 



-a. 



(A4) 



The very existence of the Thomas and Efimov effects 
clearly reveals that cold collisions of three Bose atoms are 
sensitive also to the spatial range of the interaction. At large 
distances r the asymptotic form of the inter-atomic potential 
is given by V(r) ~ -C^/r 6 , where is the van der Waals 

r— *oo 

dispersion coefficient. We shall thus determine the effective 
binary potential of alkali atoms in such a way that it provides 
a straightforward access to the full T matrix and, at the same 
time, recovers the exact scattering length a as well as those 
low energy physical observables that are sensitive also to C^, 
such as, for instance, the the energies of the trapped atom pairs 
in Fig.n As shown in Ref. Ell , the C(, coefficient enters these 
measurable quantities in terms of the mean scattering length 
of Eq. ©. In our applications to energy spectra in the vicin- 
ity of a zero energy resonance the mean scattering length a 
determines the binding energy of the highest excited diatomic 
vibrational state (cf. Fig.^ at positive scattering lengths a by 
Eq. 0. 



d. Solution of the Lippmann-Schwinger equation 

To efficiently solve the two-body and three-body 
Schrodinger equations, it is convenient to choose a (non- 
local) separable potential of the general form (see, e.g., 

Ref. stum) 



Vsep = \g)A{g\ 



(A5) 



formally, by iteration, in terms of its Born series: 



tep(z) = J] [VsepGjrW]' V sep . (A7) 

;=o 



A simple derivation then shows that the T matrix f sep (z) asso- 
ciated with y sep is given by the formula: 



T sep (z) = |g>Tftee(z)<g|. 



(A8) 



Here the function Tf ree (z) can be determined from a geometric 
series to be: 



Tfreefe) = [A' 1 - (g\G f ™(z)\g)\ 



e. Adjustment of the separable potential 



(A9) 



Evaluated at zero energy, Tf ree (z = 0) is related to the s wave 
binary scattering length a through Eq. JA4V Given the Gaus- 
sian form of \g) in Eq. iA6\ . a spectral decomposition of the 
Green's function G[j ee (z) in terms of plane wave momentum 
states shows that its matrix element in Eq. iA9\ can be eval- 
uated at zero energy to be (g|Gg lee (0)|g) = -2mcr 2 /fi 2 . This 
leads to the relation 



A = 



-H 2 /(2mcr 2 ) 
1 - y/ncr/a 



(A10) 



which can be used to eliminate the amplitude A in favour of 
the range parameter cr and the scattering length a. 

The single pole of Tf lee (z) at positive scattering lengths de- 
termines the energy E\, of the highest excited near resonant 
vibrational bound state (cf. Fig. [Q. The adjustment of the 
energy E\, to Eq. Q has been performed in Ref. ll39il and de- 
termines the remaining unknown range parameter to be: 



cr - yfna/2. 



(All) 



In our application to 85 Rb we use cr = 69.58 a.u. which cor- 
responds to C(, = 4703 a.u. ll45ll . The amplitude A depends 
on the magnetic field strength B via the scattering length a in 
accordance with Eq. ( IA101 



as an effective replacement of the full microscopic binary in- 
teraction V(r). Here \g) is usually referred to as the form factor 
that sets the scale of the spatial range of the potential, while 
the amplitude A determines the interaction strength. For con- 
venience, we choose the form factor to be a Gaussian function 
in momentum space |39|: 



^ ) = ( p|, ) = ._ 



exp 



2 1 
p cr 



(A6) 



To adjust the amplitude A and the range parameter cr in such 
a way that V sep recovers the scattering length a as well as the 
mean scattering length a of the microscopic interaction po- 
tential V(r), we determine the full two-body energy spectrum 
associated with the separable potential via its transition ma- 
trix. We thus solve the Lippmann-Schwinger equation ( IA3I 



2. Energy levels of a trapped atom pair 

a. Energy levels of the relative motion in the absence of an 
inter-atomic interaction 

In the following applications we describe the micro-trap by 
a three dimensional spherically symmetric harmonic poten- 
tial. The linear confining force then allows us to separate the 
centre of mass motion from the relative motion of a trapped 
atom pair I24ll25ll26ll . In the absence of an inter-atomic inter- 
action the Hamiltonian of the relative motion is thus given by: 



Ho 



2(f) 



V + 



1 lm\ 
2\2> 



2 2 



(A12) 
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Here a>h is the angular trap frequency. Throughout this ap- 
pendix we choose energy states \<fki m ,) of the harmonic os- 
cillator with a definite orbital angular momentum, where I is 
the angular momentum quantum number and m; is the orien- 
tation quantum number. The associated energies are given by 
Eki = hti>ho(2k + 1 + 3 12), where k = 0, 1,2,... labels the vibra- 
tional excitation of the atom pair 16311 . We denote the spheri- 
cally symmetric vibrational energy states by \ipk) - \ fkoo) and 
their energies by = E^o- Their wave functions are given by 



Eq. ( IA16I > is given by: 



<Pk(r) 



Ai) 



where L k 2 is an associated Laguerre polynomial. The param- 
eter f3 = mu)h l{2h) is related to the harmonic oscillator length 
«ho for a single atom, i.e. the trap length, by f3 = l/(2a? ). 



b. Separable potential approach in the presence of a trapping 
potential 

The spherical symmetry of the trap allows us to determine 
the energy levels of the relative motion of a pair of interacting 
trapped atoms in complete analogy to their counterparts in free 
space The associated Hamiltonian is then given 

by 



(g\G (z)\g) = J] 



k=0 



l<Wlg)f 

z-E k 



(A17) 



During the course of our studies, we have compared the sep- 
arable potential approach to the two-body energy spectrum to 
predictions obtained with a microscopic potential V(r) for a 
variety of scattering lengths and trap frequencies. Figure ^ 
shows such a comparison for a rather tight Vh = 300 kHz 
atom trap, which clearly reveals the applicability of the sep- 
arable potential approach in the range of energies relevant to 
the adiabatic association of Efimov trimer molecules. 



APPENDIX B: THREE-BODY ENERGY LEVELS AND 
WAVE FUNCTIONS 



In this appendix we derive the Faddeev equations that deter- 
mine the exact energy levels of three interacting Bose atoms 
in the confining potential of a spherical trap. We then describe 
a practical method to exactly solve these equations in the sep- 
arable potential approach. 



H = H + V(r), 



(A14) 



1. Faddeev approach 



where V(r) is the spherically symmetric microscopic inter- 
atomic potential and Ho is the Hamiltonian of Eq. (IA12I . In 
the following, we shall denote the energies associated with H 
by E 2 b- 

Given that typically the trap length Oho very much exceeds 
the van der Waals length Z V dW> the full microscopic interaction 
can be replaced by the separable potential of Appendix IA ll to 
describe the limited range of energies involved in the adiabatic 
association of molecules. The associated T matrix T sep (z) of 
the relative motion of a pair of trapped interacting atoms can 
then be determined in analogy to Appendix IA II This yields: 



\g)T(z)(g\. 



(A15) 



The function r(z) can be obtained from Eq. ( IA9> by replacing 
the Green's function G Q Ke (z) in free space by its counterpart 
Go(z) = (z - Ho) 1 in the presence of the trapping potential. 
This leads to the relation: 



r(z) = 



(TfreeWr 1 " 4 CT 2 /h 2 - (g\6 (z)\g) 



(A16) 



The poles of r(z) determine the energy levels of the Hamilto- 
nian JA14> in the separable potential approximation, i.e. the 
poles are located at the energies z = Erg- 

We have evaluated the function r(z) in terms of the oscilla- 
tor states \<pk) of Eq. jA13> using the spectral decomposition 
of the Green's function Go(z). The matrix element relevant to 



a. Three-body Hamiltonian and Jacobi coordinates 

Throughout this appendix, we assume that the atoms inter- 
act pairwise via the potential V(r). The complete Hamiltonian 
is then given by: 



H = H + V(r 23 ) + V(r 3 i) + V(r 12 ). 



Here the non-interacting Hamiltonian 



(Bl) 



Ho 



Zl + -muZx] 



(B2) 



accounts for the kinetic energy and the harmonic trapping po- 
tential of each atom, while V(r,j) describes the interaction be- 
tween the atoms ; and j = 1, 2, 3) in dependence on their 
relative coordinates ry = r, - tj. In the following, we employ 
the Jacobi coordinates R = (ri +r2 + r3)/3,p = ri-(r2 + r3)/2 
and r = r 2 - r 3 of Fig.|S]to separate out the three-body centre 
of mass. In analogy to the case of a pair of trapped atoms, the 
harmonic force then allows us to divide the non-interacting 
three-body Hamiltonian of Eq. JB2t into three harmonic os- 
cillator contributions associated with the Jacobi coordinates. 
The binary potentials involve only the relative coordinates p 
and r. The complete three-body Hamiltonian can thus be rep- 
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resented by: 




FIG. 6: Jacobi coordinates of the relative motion of three atoms. The 
set of coordinates p and r is selected in such a way that it is suited 
to describe the hypothetical situation of an interacting pair of atoms 
(2, 3) with atom 1 playing the role of a spectator. 



b. Faddeev equations for three trapped atoms 



the three-body energy state is given by = + \i//2) + 1*/^)- 
Inserting this Faddeev decomposition into Eq. ( IB5i on the left 
hand side and rearranging the terms yields: 

[1 - G (£)Vi] |<A,> = G (£)Vi (|*A 2 > + l<A3» . (B7) 

Equation (IB7> can then be solved formally for \if/\) by mul- 
tiplying both sides with [1 - Go(£')Vi] _1 from the left. This 
leads to the Faddeev equation: 

|<Ai> = [1 - GoiEWiV 1 G (E)Vi (\iff 2 ) + l<A 3 » • (B8) 

The kernel of Eq. (IB 8 i can then be expanded into the power 
series associated with the inverse matrix [1 - Go(E)Vi]~ l . 
This expansion yields: 

00 

[1 - GoCEmr 1 G (E)Vi = Go(£)J] [ViG (E)V V x . (B9) 

7=0 

In analogy to Eq. JA7L the sum on the right hand side of 
Eq. dB9l can be identified as the Born series associated with 
the Lippmann-Schwinger equation: 

Ti(z) = Vi + V x G Q (z)T x (z). (BIO) 

Representing Eq. dB8b in terms of the T matrix T\(E) for 
the interacting pair of atoms (2, 3) then recovers the Faddeev 
equation for in its original form 13111 : 

\i/M) = G (E)T 1 (E)(\i/, 2 ) + \i/ fi )). (Bll) 

The Faddeev equations for I1//2) and 1 1^3) are obtained by cyclic 
permutations of the atomic indices. The resulting set of three 
coupled matrix equations for an d I1A3) determines all 

energy levels of three interacting atoms in a trap. We note that 
the formal derivations leading to Eq. dB 1 1> do not refer to the 
specific nature of the confining potential. This is the reason 
for the general validity of the Faddeev approach in free space 
as well as in the presence of an atom trap. 



The Faddeev equations for the energy levels of three pair- 
wise interacting atoms in the confining potential of a trap 
can be derived largely in analogy to their counterparts in free 
space 13111 . To this end, we introduce the Green's function 

G (z) = {z- HoT 1 (B4) 

associated with the non-interacting Hamiltonian of Eq. dB2t 
and denote the binary potential associated with the atom pair 
(2, j) by Vk — V(Tij) for each of the three possible combina- 
tions of atomic indices 

(i/t) = (l,2,3),(2,3,l),(3,l,2). 

The stationary Schrodinger equation H\ip) = E\ifS) for a three- 
body energy state \if/) can then be represented in terms of the 
matrix equation: 

\i f ) = G (E)(V l + V 2 + V 3 )\if). (B5) 

Introducing the Faddeev components 

m = G (E)Vi\iff), (B6) 



c. Faddeev approach for three identical Bose atoms 

In the special case of systems of three identical Bose atoms, 
such as S5 Rb3, the three Faddeev components depend on 
each other through cyclic permutations of the atoms. These 
permutations can be represented by a unitary operator U, 
i.e. \\j/2) = U\if/i) and I1//3) = 1/ 2 |i//i), which transforms the 
three-body wave functions in accordance with the formula: 

('W^)(R > p > r) = ^(R 1 p',r'). (B12) 

Here the primed coordinates are determined by {^,^ = JTl,), 
where 

(B13) 

is a 6 x 6 matrix satisfying JT 3 = 1 ■ The first Faddeev com- 
ponent of Eq. JB 1 ll > is thus determined by the single Faddeev 
equation: 

|<Ai> = G {E)Ti{E) + U 2 ) l<Ai). (B14) 
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d. Basis set expansion approach 

In the following, we employ a basis set expansion approach 
to solve the Faddeev equation JB14> . which provides an exten- 
sion of the momentum space Faddeev approach 15^1 to sys- 
tems of trapped atoms. Since, according to Eq. JB3i . the free 
Hamiltonian Hq can be divided into a sum of three indepen- 
dent harmonic oscillators, we choose the basis of products 

\^>KLM L ,lpKAiA A ,iPklmi) = \^KLM L )\t/>KAfi A )\<Pklm,) (B15) 

of the energy states of the individual oscillators with a definite 
angular momentum, i.e. \<& K lm l ), and \<p Umi ) are the 

oscillator energy states associated with the Jacobi coordinates 
R, p and r, respectively. 

The kernel of the Faddeev equation JB14b is diagonal in the 
non-interacting energy states \$>klm l )- The component \if/\) 
can thus be chosen in such a way that it factorises into a centre 
of mass part and a relative part, i.e. 

I-Ai) = VbKuiJWf). (B16) 

The relative part li/'" 1 ) then satisfies a reduced Faddeev equa- 
tion at the shifted energy Ej, B = E — Ekl- 

|iAf > = Gf (£ 3B ) Tf (£ 3B ) \U + I/ 2 ] |^ eI >. (B17) 

Here Gjj el (£36 ) and T[ el (Et,b ) can be interpreted in terms of 
a reduced non-interacting Green's function and a reduced T 
matrix, respectively, which depend only on the Jacobi coordi- 
nates p and r describing the relative motion of the three atoms. 

In order to solve Eq. JB 17i . the basis set expansion ap- 
proach takes advantage of the convenient diagonal representa- 
tion of the reduced non-interacting Green's function in terms 
of the chosen basis states: 

G reI CE 3B ) V V V ^ KA ^ ,tfklm '^ KA ^ , ' fklm ^ 
° 3B ~ JlfoAA E, B -E KA -E kl 

(B18) 

Since the reduced T matrix rj el (£33) includes just the inter- 
action between the pair of atoms (2, 3), it is related to the T 
matrix T(z) of the relative motion of this atom pair by the for- 
mula: 

(4 , Kfy i ,(Pklmi\Tf l (E3B)\4>KU'ii' ;l ,<Pk'l'm' l ) =(<Pklmi\T(z)\ipk'l'm' l ) 

Here z = Ei, B - E kA accounts for the energy of atom 1 . The 
kernel of Eq. JB14I > is thus completely determined by the so- 
lution of the two-body Lippmann-Schwinger equation in the 
presence of the trapping potential. 

The complete three-body energy state can be factorised in 
analogy to Eq. JB 16i . i.e. 

W) = \<&KLM L )\lp K \)- (B20) 

Given the solution of Eq. JB17> . the reduced Faddeev compo- 
nent \4> r f) determines |^ re i) by the relationship: 

| i A re i> = (l+1/ + 'Z/ 2 )|<>. (B21) 



Equations JB 17i . JB 1 8i and JB19i set up our general approach 
to the energy spectrum of three interacting atoms in a trap, 
while Eqs. JB20> and (IB2U yield the associated three-body 
energy states. 



2. Solution of the Faddeev equations in the separable potential 
approach 

a. Faddeev equations in the separable potential approach 

We shall show in the following that the separable poten- 
tial approach to the two-body T matrix, in combination with 
the basis set expansion, provides a practical scheme to exactly 
solve the Faddeev equation JB17> in the presence of a trap- 
ping potential. To this end, we apply Eqs. JB 19i and JA15> to 
EqjBTT}. This yields: 

CO CO A 

l^ 1 ) =Gf (E 3B ) J] Z Z \^8)r(E 3B - E KA ) 

k=0 A=0 fi*=-A 

x <0,^„^| + I/ 2 ] (B22) 
Equation JB22I > reveals that li/^ 6 ') is of the general form: 

\ift 1 ) = G% l (E 3B )\f,g). (B23) 

As only the unknown amplitude f K ^ /lA = needs to 

be determined, the separable potential approach significantly 
simplifies the Faddeev equation JB 17i by reducing its dimen- 
sionality from six to three. We shall show in the following that 
due to the spherical symmetry of the form factor \g) the num- 
ber of dimensions reduces even further to only one provided 
that the three-body energy levels under consideration have s 
wave symmetry. 

b. Basis set expansion 

In accordance with the basis set expansion approach, 
we consider the projected amplitude f K Aft r The ansatz of 
Eq. JB23> for the solution of the Faddeev equation JB22I > then 
determines the amplitude f K x fli by the matrix equation: 

OO DO A' 

fnAUt = T(Ej B - E kA ) ^ ^ ^ I KkA/i a ,k'A'ii' 1 (E3b) fK'A'/i' A - 
K'=0A'=0fi\=-A' 

(B24) 

In accordance with Eq. JB 181 . the reduced kernel matrix asso- 
ciated with this equation for f K/ ^ A is given by the formula: 

%awfX4 (E3b) = (4> k a^,8\ (U + U 1 ) Gf (E 3B ) \cf> K , A% ,g). 

(B25) 

The complete kernel also involves the function t(Ej, b - E K/X ) 
which we have discussed in detail in Appendix IA II Inserting 
the spectral decomposition of the Green's function of three 
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non-interacting trapped atoms of Eq. JB 1 8i then determines 
the reduced kernel matrix to be: 

k,k<=0 U'=0 »„=-/,«;=-/' 3B C ^ £ ' A / 
X ((pKfy^tpklmtl \U + I/ 2 ) {(pK'A'tt'^fk'l'm^- 

(B26) 



c. Symmetry considerations 

The spherical symmetry of the form factor |g) of Eq. iA6\ 
implies that only the spherically symmetric basis states \<pk) = 
\fk(a) contribute to the summation in Eq. ( IB26> . i.e. 1 = 1' = 
and mi = m', = 0. As, moreover, the total angular momentum 
operator X = A + 1 associated with the three-body state \4> r e\) 
commutes with the permutation operator K, the off-diagonal 
elements X + X and p. + // of the kernel vanish, i.e. the 
kernel does not couple solutions of different angular momenta 
X. Furthermore, the Hamiltonian H commutes with 11 which 
implies that the matrix element of {^U + K 2 ^ in the reduced 
kernel matrix of Eq. ( IB26> is non-zero only if 



k + k = at' + k'. 



(B27) 



This is a consequence of energy conservation. 

We shall focus in the following on those energy states of 
three trapped interacting atoms that correlate adiabatically, in 
the limit of zero trap frequency, with the three-body s wave 
Efimov states. We thus restrict the discussion to three-body 
states with zero total angular momentum £ = 0. In this case 
all angular momentum quantum numbers are zero and may 
be omitted. This implies that the kernel matrix in Eq. {524} 
reduces to 

A gki^vriiu + f U 2 )\4>^y k >)gk' 
7Ck'0E3b)= 2j ^ > ( B28 > 



k,k'=0 



E3E 



and the amplitude f K = {(p K \f) satisfies the matrix equation: 

00 

f K = t (E 3B - E K ) J] <K KK , (E 3B )f„ . (B29) 



K'=0 



In accordance with the ansatz of Eq. JB23> . the first Faddeev 
component is then given in terms of the amplitude f K and the 
form factor gk by the formula: 



f«gk 



E3B - E K - Ek 



(B30) 



3. Determination of the kernel matrix 

A main difficulty in the numerical determination of the 
amplitude f K from Eq. ( IB29t consists in calculating the re- 
duced kernel matrix K^* (£33) for the variety of trap fre- 
quencies studied in this paper. While in tight atom traps 



(vho > 100 kHz) the discrete nature of the energy levels is 
most significant, the opposite regime of low trap frequencies 
(vho < 1 kHz) involves a large range of vibrational quan- 
tum numbers k leading to a continuum of modes in the limit 
Vho — * 0. To obtain a stable scheme for the determination 
of the reduced kernel matrix in the limits of both high and 
low trap frequencies, we have performed separate treatments 
of the regimes of low and high vibrational quantum numbers 
k. Since these aspects of the studies of trapped systems of 
three interacting atoms differ significantly from the known 
techniques to solve the Faddeev equations in free space l56ll . 
we shall outline in detail the numerical procedure we have ap- 
plied. 

a. The kernel matrix at low vibrational excitations 

In the limit of small k the determination of the reduced ker- 
nel matrix < K KK i consists mainly of the calculation of the matrix 
elements 

C Kk yk> = (<t>«, fk\ \U + <U 2 \ \<t>«,tp k -)- (B31) 

We perform this calculation in the configuration space repre- 
sentation. To this end, it is convenient to introduce the auxil- 
iary function: 



\j 2n T(k + §) 

This function is related to the spherically symmetric harmonic 
oscillator states (cf. Appendix IA 2> associated with the Jacobi 
coordinates p and r by the formulae 



<fik(r) =i/,k(J3 r ;P r r 2 ), 



(B33) 
(B34) 



respectively. Here the parameters f3 p = |ma»h /^ and f3 r = 
jmoJho/fr account for the different masses associated with the 
individual harmonic oscillator contributions to the three-body 
Hamiltonian ( IB3i . The matrix element involving the permu- 
tation operator can then be represented by: 

C*xt!= j d 3 pd 3 ri(r K (/3 p ;/3 p p 2 )^k(/3 r ;l3 r r 2 ) 

X [tf /K ,(J3 p -/3 p p' 2 )tfrk > (J3 r -/3 r r' 2 ) 

+ ^(Pprfpp'^k-iPrlPrr" 2 )]. (B35) 

Here the permutation operators K and H 2 transform the co- 
ordinates p and r into primed and double primed coordinates 
= J\r) an d (r") = ^ 2 (r)' respectively, in accordance with 
the transformation matrix JT of Eq. JB 1 31 . The transformation 
with J' yields: 



PpP' 2 =\p P P 2 + \p r r 2 - fypPr P ■ r, 



JlJa .2 , 1 2 /3. 



0rT" = ^ P p Z + -fiS + y -p p pr P ■ r. 



(B36) 
(B37) 
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This implies the relationship: 

P P p' 2 +p r r' 2 =/3 P p 2 + p r r. 



(B38) 



Similar relations hold for the double primed coordinates with 
the reversed signs in front of the square roots. The integrand 
on the right hand side of Eq. dB35b . therefore, depends only on 
the variables p and r in addition to the variable x — p ■ r/(pr), 
which involves the angle between the coordinates p and r. 
The integration over the remaining three variables is readily 
performed and leads to the formula: 



p 2 dp 



r 2 dr ^ K (J3 p -p p p 2 )^ k {p r \P r r 2 ) 



i: 



dx i/v 



3 ~ .2 , 1 2 



PA ~P p p L + -py + yj-P p P r prx 



PpPrprx 

J 

(B39) 



This formula can be further evaluated by changing the vari- 
ables to u = p p p 2 and v = p r r 2 and using the explicit form 
of the harmonic oscillator wave functions of Eq. (IB32I . This 
evaluation yields: 



/ r(x + i)r(fc + pr(x' + i)r(*' + i) 



r(K+hT(k+hr(K' + hT(k' + h 



(B40) 



The coefficients I KK ,K'k' that involve the associated Laguerre 
polynomials read: 



UK > k , = I du uie~"L K 2 (u) I 
Jo Jo 



dvvi e - v l}}\v) 



x: 



dx P K 'k'(u, v, x). 



(B41) 



Here the function 



P K . k .(u,v,x) =L 



1 3 

— It H V — 

4 4 



3 

—uvx 
4 



rii) 3 



1 



x LJ \-u + -v -. 
* 4 4 V4 



3 

—uvx 



(B42) 



is a polynomial of degree K'+k' in the variable x. These deriva- 
tions reveal that the matrix element in Eq. JB40> is indepen- 
dent not only of the inter-atomic interaction potential but also 
of the frequency of the atom trap. The trap frequency enters 
Eq. ( IB29> through the projections of the form factor onto the 
basis states and through the energy denominator. 

To further evaluate Eq. JB4H . we represent the function 
P K 'k'{u, v, x) by the sum: 



K'+k' 



PMu,v,x)=YjPlk>( u > v ) xy 

7=0 



(B43) 



over the variable x in Eq. (IB41i . Only the even powers x y of 
the variable x contribute to this integral, while all terms in- 
volving the odd powers vanish. When y is even, however, it 
turns out that the coefficients P Y K , k ,(u, v) themselves are bivari- 



ate polynomials in the variables u and v and can also be ex- 
panded in powers of u and v 
This expansion thus yields: 



panded in powers of u and v with expansion coefficients P y J' k , . 



K'+k' K'+k'- 



(B44) 



T=0 /=() 



The representation of the expansion coefficients P r K , k ,(u,v) 
in Eq. JB44> in terms of a polynomial allows us to take advan- 
tage of the general formula 



f 

Jo 



(i) I : < t < k 

= \ (_i)<T( T + |)( T )- k< t 



(B45) 



to perform the remaining integrations in Eq. dB41i over the 
variables u and v. Here Pj is a combinatorial. Equation JB45l i 
thus determines the coefficients I^'V of Eq. JB4U to be: 



K' + k' r. K'+k' J -, \ I \ 

y=Q * T =K V f \ / 

"WOP- 



(B46) 



The limits of the different sums in Eq. ( IB46> lead to a fur- 
ther simplification in the determination of I K k, K 'k' as follows: 
The summation over the index t is limited by the condition 
k' + k' - t > k, which, together with the condition r > k for the 
index r, implies the inequality K'+k' > K+k. Similar consider- 
ations show that also the inequality K'+k' < k + k is fulfilled, 
which, in summary, leads to the restriction K + k = k' + k' . 
These explicit derivations simply recover Eq. JB27L which 
we have obtained independently from general symmetry con- 
siderations. Consequently, the summations over t and t in 
Eq. JB46> reduce to a single term determined by r = K and 
t — k. This gives the coefficient I KK yk' to be: 

W* = (-ir*r(ir + |)r(* + |) 't^Kt CB47) 

(y even) 

The remaining, as yet, undetermined expansion coefficients 
P 7 J k , of Eq. (IB44> can be obtained from the explicit form of the 

associated Laguerre polynomials, i.e. L K 2 (x) = 2y=o c kvX v , 
and their expansion coefficients: 



Ckv = (-D V 



1 l K + h 



(B48) 



Here the coefficients P 7 .,,(u, v) depend on the variables u and 
v. Equation (IB43I then allows us to perform the integration To this end, we consider the function P y . k ,(u,v) of Eq. (IB43> . 
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whose explicit form can be determined from Eq. JB42> to be: 

k' min(v,y) 



pi k xu,v) ( M v)«2> v j>„ £ ( -ir 

^ v=0 ti=0 a=max(0,y->i) 



4" + 4 vl 



(B49) 

Equations (IB44> and ( IB49> then give the remaining expansion 



coefficients P^, to be: 



p!' 



min(/c',y) 



=4^ " k C K , K ,Ck>k> 



a=max(Q,y-k') 
min(A:'— o-,A"-^) 

z 



k' 

y — a 



ar\ / k f + a - y 



^=max(0,/f / — /f— cf + ^ ) 



(B50) 



In order to summarise our results in a concise form, we 
introduce the coefficients 



pjKk 
r K'k' 



p yKk 
r K'k' 

C k' K'Ck'k' 



(- 1 y +k rv + i)r(&' + 1 )p y t , (B5 1 ) 



which are obtained simply by dividing P y *t of Eq. JB50> by 
the expansion coefficients c K i K > = (-1)* /HY + 1) and = 
(-1)*' /r(k' + 1) of the associated Laguerre polynomials. An 
analysis of its summation limits shows that the sum over 77 in 
Eq. ( IB50> vanishes unless the conditions y < 2k and y < 2k 
are fulfilled. In accordance with Eq. (IB47> . we thus obtain: 



r(/f+f)r(fc+§) ™n(K+u*,2*) 
hk,K'k' = TTJ-TT TTT^TTT - TT 2-t 



T( K ' + i)r(k' + 1) 



7=0 

(7 even) 



7+1 



(B52) 



The complete matrix element of Eq. JB35b is then determined 
by the formula: 



rp+i) Htt + f) r(k+i) r(fc + f) 
~\ t{k' + 1) r(^' + h iXJf + 1) r(jf + 1) 



min(/H-&,2/f,2£) 



S 

7=0 

(7 even) 



2 - 



7+1 



ykk 
K'k'' 



(B53) 



Since the triple summation over a, r\ and 7 in Eqs. JB50> 
and JB53i contains only integer numbers, the calculation of 
the matrix element C^^'k' is, in principle, straightforward. As 
these numbers can, however, be very large in magnitude and 
have alternating signs, we have employed a computer algebra 
system to carry out the sum in multi-precision integer arith- 
metic. This proved practical for the indices < k, k' < 40. 



b. The kernel matrix in the limit of high vibrational excitations 

The numerical determination of the matrix element C K kyk' 
using Eq. dB53l > becomes impractical in the limit of large in- 
dices k and k' . We shall, therefore, provide a scheme to di- 
rectly determine the reduced kernel matrix for these higher 



vibrational quantum numbers in terms of its asymptotic form 
in the continuum limit. Starting from Eq. ( IB25> evaluated at 
zero angular momenta A — and X = 0, the exact reduced 
kernel matrix is given by the formula: 



?C«< (E 3 b) =(<?>k, g\ (14 + U 2 ) (£ 3B ) \<f>K',g) 

{<l> K ,g\in4 + f U 2 )\<i> K ',<f>k)(<f>k\g) 



=1 

k=0 



E3B ~ E K > - E K 



(B54) 



It turns out that, due to the permutation operators, the right 
hand side of Eq. JB541 depends just on the oscillator wave 
functions <p K (fi), <t>K'{p) and (f>k{r) in a limited range of radii 
p and r set by the width of the form factor g(r). The typi- 
cal length scale associated with the width of g(r) is set by the 
range parameter cr of Eq. JA1 1> and is thus determined by 
the van der Waals length Z v dw- As the van der Waals length 
is typically much smaller than the harmonic oscillator length 
flho = VS/( mw hoX the trapping potential is flat within the 
relevant range of radii p,r < cr and the potential energy is 
much smaller than the mean kinetic energy of the highly ex- 
cited oscillator states. We can, therefore, replace in Eq. (IB54> 
the Green's function G\f (E^) of the trapping potential by its 
counterpart G[j ee (£36 ) in free space and perform the contin- 
uum limit of the oscillator states \<p K ) and \<p K ')- 

Given that the radius p is limited by the condition p < cr, 
in the continuum limit the harmonic oscillator energy wave 
functions <f>K,i^ A (p) approach, up to a normalisation constant, 
the partial waves \qA/i,\), i.e. the improper energy states of the 
free space Hamiltonian 



h 2 



2(2m) 



(B55) 



with a definite angular momentum. 

These partial waves thus satisfy the Schrodinger equation 
H p \qA/u,i) = E q \qAnx), associated with the kinetic energy E q = 
# 2 /[2(|m)], and are related to the plane waves |q') by 



<q'|?4u> = (~i) A qMq ~ <?')C (q'/V), 



(B56) 



where T^ 1 is a spherical harmonic. 

To determine the asymptotic form of the reduced kernel ma- 
trix in the limit of high vibrational excitations, we insert into 
Eq. JB54> two complete sets of improper states \qA^i^) and 
\q' A'ii'j). This yields: 

KkAEw) ~ i dqq 2 (cf> K \q00) I dq' q' 2 {q'00\(f>K') 

K,K'^CO J J 

x faOO, g\ (<U + 1/ 2 ) G ree (E 3B ) \q'00, g). 

(B57) 

A simple calculation then shows that for high vibrational 
quantum numbers k the function q 2 ((p K \q00) is sharply peaked 
about the central momentum q at the matched energies E K = 
E q . This implies: 



4 * L 3 
q - \l —mncoix, \2k+ — 



(B58) 



17 



A similar relation holds for the central momentum q' 
associated with the function q' 2 (q'00\(f> K i). We may, 
therefore, evaluate the slowly varying matrix element 
{q00,g\ (t{ + W 2 )G f ™ (E m )\q'00,g) in Eq. (|B53 at q = q 
and q' — q' . The remaining integration over q then yields 



by the formula: 



I 

Jo 



dqq 2 (^\q00) 



t2mhuh '\' 1 ' 

3 



r(*+§) 
t(k+\y 



(B59) 



and the integration over q' can be performed in an anal- 
ogous way. Using the spectral decomposition of the free 
space Green's function in terms of plane wave momen- 
tum states and the explicit Gaussian expression for the 
form factor in Eq. (IA6> . the remaining matrix element 
{qOO, g\ (<Z/ + %l 2 ) Gf[ ee (£ 3 b) l^'OO, g) can be determined an- 
alytically. In the limit of high vibrational excitations k and k' 
the reduced kernel matrix is then given by the formula: 



„ , 16 |r(^+|)r(^' + |)/2m Who cr 2 



^ \ro + i) T(k' + i) \3 n 



,2\3/2 



m 

x — — exp 

qq' 



-(q 2 + q' 2 ) - mE 3B 



cr 2 /h 2 



Ei [(mE m - (q 2 + q' 2 - qq'))o- 2 /h 2 ] 



- Ei [(mEm - (q 2 + q' 2 + qq')) o- 2 /h 2 j 

(B60) 

Here Ei denotes the exponential integral Ei(x) = dt r l e'. 

4. Three-body energy wave functions 

a. Basis set expansion of a three-body energy state 

Once the reduced kernel matrix < K KK '(Et >b ) has been calcu- 
lated, the three-body energies E^ and amplitudes f K = {<p K \f) 
can be obtained from the numerical solution of Eq. ( IB29I . Ac- 
cording to Eqs. JB21I and ( IB23I >. each solution completely de- 
termines its associated three-body energy state by the formula: 



\cf> K ,<p k l£=0,Mj:=OJJ)= j] 



(-D 



I-)!!/ 



i, ^^2^+~l 



\4>Khn,,<Pkl-mi)- 



(B62) 

Here (-1)'-""/ V27TT is a Clebsch-Gordan coefficient. Al- 
though the permutation operator 14 commutes with the to- 
tal three-body angular momentum _£, it does not individu- 
ally commute with the partial angular momenta A and 1. De- 
spite the fact that \ f) and \g) contain only contributions from 
the spherically symmetric harmonic oscillator basis states \<p K ) 
and \ipk), respectively, the basis set expansion of the complete 
state |tArei) thus involves oscillator states \<p K A^ A ) and \ipki mi ) with 
all A, ^ = 0, 1, 2, . . . values. This leads to the difficulty that 
the matrix elements of Eq. dB35l need to be calculated also 
between these states. While this effort will generally be in- 
evitable, we shall show that, for the purposes of this paper, 
it can be avoided. We utilise the fact that the determination 
of matrix elements, in the three-body energy states, of those 
quantities that commute with all three operators 14, A, and 1 in- 
volve just the spherically symmetric oscillator states \<p K ) and 
\(fk). This is most evident for the normalisation constant and 
the orthogonality relation of the three-body energy states. 



b. Orthogonality and normalisation 

To derive the normalisation constant of the fully sym- 
metrised three-body energy states as well as their orthogonal- 
ity relation, we consider a pair of three-body states |^ re i) and 
l^rei) with the associated energies £33 and £33, respectively. 
Starting from Eq. (IB6H and using the identity 14 3 = 1, the 
overlap between these states is given, in terms of their first 
Faddeev components, by the matrix element: 

<^re# re i> = 3<<'| (1+04 + 14 2 ) liAf). (B63) 



Since both \if/ r f) and \if/ r f) can be expanded in terms of the 
spherically symmetric oscillator states \<p K ) and \<pk), we obtain 
the formula: 



|^ re i> = (l+1/ + 1/ 2 )G^(^3B)l/,g>. (B61) .— V/Trelu \ 



The numerical determination of the complete state \>p m \) re- 
quires an expansion of Eq. JB6H into harmonic oscillator 
states. Since, throughout this paper, we focus on three-body s 
wave states, a suitable basis set for the expansion is provided 
by the the three-body oscillator energy states with zero total 
angular momentum £ — A + 1. The total angular momentum 
quantum number of a three-body s wave state is thus given 
by Si — 0, and Aix = is its associated orientation quantum 
number. The s wave basis states are then related to the oscil- 
lator energy states and \ifki m ,) with A — I and //,) = -mi 



:,k=0 
K+k 



{(j) K ,iPkWi x ) 



(B64) 



K'=0 



Here we have used Eq. ( IB27> to eliminate the summation over 
k! in favour of the relationship k' = k + k — k ' . Inserting 
the ansatz of Eq. (IB23> to eliminate the Faddeev components 
[^ re i) and |i^ re i) in Eq. (IB64> in favour of |/) and |/), respec- 
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tively, then leads to the orthogonality relation: 

f*gk 



<^re#rel> =3 ^ — 



^r (£3B -E K - E k )(E jB -E K - E k ) 

K+k 

f«gk + Yj^" + <Ut ) 1^' <Pu)f*8K 



k'=0 



(B65) 



In the special case of |i^ re i) = ItArei) Eq. JB65i also yields the 
normalisation constant of the three-body energy state \if/ K \). 



c. Hyper-radial probability density 

A symmetrised zero angular momentum three-body wave 
function i// K \(p,r) = (p, r|^ re i) depends only on three vari- 
ables among the six Jacobi coordinates @|. These variables 
may be chosen, for example^as the radii p = \p\, r = |r|, and 
the angle between p and r 16411 . A common way to visualise 
the spatial extent of three-body states by a one dimensional 
function involves the transformation of the Jacobi coordinates 
p and r to hyper-spherical coordinates. Among these coordi- 
nates the hyper-radius and the hyper-angle are given by 



[m 2 „ 1 „ 

V yUR V 3 2 



(B66) 



and tan <D(p, r) = y| - f , respectively. Apart from R all hyper- 
spherical coordinates are angular variables. The "mass" pa- 
rameter fiR in Eq. JB66> ensures that R has the unit of a length. 
In Fig. [5] we have chosen it to be fiR = m. The hyper-radial 
probability density associated with the three-body state |^ re i) 
is determined, in terms of the projection operator 



d 3 pd 3 r\p,r)6(R-R(p,r)) (p,r\, 



by the formula: 



P{R) = (MPrIM- 



(B67) 



(B68) 



For a normalised three-body state P(R) satisfies the normali- 
sation condition L dR P(R) = 1 and provides a measure of 
the spatial extent of its wave function. 

As the hyper-radius of Eq. JB66> does not depend on the 
angular variables and is also invariant with respect to H, 
the probability density P(R) can be calculated, similarly to 
Eq. JB64> . just in terms of the spherically symmetric oscilla- 
tor states \<p K ) and \tpk). This yields 



p(R) =3 Y^i^PR^n) 



K+k 



2<0*. ip k \ [v + <u 2 ) \^,^)(<p K ,, n ^f) 



k'=0 



(B69) 



where k' is given by the relationship k' - K + k- k' . To deter- 
mine the matrix element 

<^V*|<fe^> = J d 3 pd 3 r [<(p,r)]* 

xd(R-R(p,r))<f> K (p)<p k (r), (B70) 



we first represent the delta function in Eq. JB67> by the for- 
mula: 



6(R 



ol/2 

- R(p, r)) = ^- 8(p; li2 {p R R 2 -P p p 2 ) 111 - r) . (B71) 

Here we have introduced the parameter /3r = p.R<x>h lfi. We 
then analytically perform the five integrations over r and over 
the solid angles associated with p and r and substitute in the 
remaining integral the variable p in accordance with /3 p p 2 = 
I3rR 2 w 2 . To represent the matrix element in Eq. JB70> in a 
concise form, we introduce the function 



f?(JiRR 2 -w 2 )= 



K\k'=0 



F(*' + 1) F(ft' + 1) f K .g k , 
I + §) F(k' + §) E 3B - E K , - E k , 

x Lfq3 R R 2 W 2 )L^ } (PrR 2 (1 - w 2 )), (B72) 

which can be evaluated numerically. Given the explicit form 
of the harmonic oscillator wave functions in Eq. jB32t , the 
matrix element of Eq. dB70b can then be obtained from the 
formula: 



{ipf\P R \<l> K , ip k ) ^P^ifiRR 1 ? 11 exp(-p^) 



/ r ( , + i)r (fc+ i) r^ wVr - 
\r(Ar + §)r(jfe+|) Jo 

x 4* WV)L<* W 2 (l - w 2 )) 
x^ l (J3 R R 2 ;w 2 ). 



w 2 w 2 



(B73) 



In this formula the integration over w lends itself to a numer- 
ical evaluation by a Gauss-Chebyshev quadrature rule of the 
second kind. 



5. Numerical implementation 

The analogue of the momentum space Faddeev approach 
ll56ll to three interacting Bose atoms in a trap involves, con- 
trary to its free-space counterpart, discrete Faddeev equations 
represented by Eq. ( IB29> in the separable potential approach. 
In the limit of low frequency trapping potentials, however, 
the basis set expansion leads to large kernel matrices and the 
numerical determination of the fixed points in Eq. dB29> be- 
comes impractical. 

In order to demonstrate the orders of magnitude of the 
kernel matrix, we consider the projection of the form factor 
gk - (<Pk\g) onto the harmonic oscillator basis. A simple cal- 
culation based on Eq. (IA6> shows that this projection is given 
by the formula: 



gk 



3/4 



n l ' 4 (i + o V2 U +(! \r(fc+i) 



(B74) 
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FIG. 7: The reduced kernel matrix "K^(E^) of Eq. <B29> . at low 
vibrational excitations k, k' < 100, for a Vho = 300 kHz trap and a 
magnetic field strength of B = 158. 1 G. The transition from the linear 
to the connecting exponential mesh at k = N e =40 is illustrated 
by the logarithmic scale of the axes associated with the vibrational 
quantum numbers k and k' (cf. Eq. <B75t ). In the linear part (0 < 
k,k' < N e ) the reduced kernel matrix was evaluated exactly using 
the results of Appendix IB 31 while the approximation of Eq. <B60> 
was used in the exponential part. The logarithmic scaling of the axes 
associated with k and k' prevents the row k = and the column k' = 
from being displayed. 
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i i i i i i 
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K 

FIG. 8: (Color online) Solutions f K of the matrix equation <B29> . us- 
ing the reduced kernel matrix shown in Fig.Q versus the vibrational 
quantum number k. The legends show the energies £33 - |/iVh of 
the associated three-body levels of interacting 85 Rb atoms relative to 
the zero point energy of hypothetically non-interacting atoms in the 
v ho = 300 kHz trap. 



Here we have introduced the dimensionless parameter £ = 
B r cr 2 = cr 2 };mu)h /ti. The form factor thus decays like gj, ~ 
e~~ lkt > in the limit of large vibrational quantum numbers k. In 
the case of 85 Rb atoms the parameter f is on the order of 
1.7 x 10~ 2 for a Vho = 300 kHz atom trap, while it gets as 
small as 1.1 x 10~ 5 at a trap frequency of Vho = 200 Hz. If we 
estimated the order of magnitude of a numerical cut-off £ max , 
for instance, by supposing that the form factor is well repre- 
sented when its value gf, at k = k mRX has decayed to 10~ 3 go 
we would require only k max as 200 basis states in the case of 
a 300 kHz trap. For a 1 kHz trap, however, this number would 
increase to k msx « 320000. At low trap frequencies the ker- 
nel matrix would, therefore, become too large for a numerical 
solution of Eq. dB29l > to be practical. 

To account for the discrete nature of the basis set expansion 
on the one hand and the exponential scaling of gk, in the limit 
of large k, on the other hand, we have introduced a linear mesh 
of length N e = 40 covering each of the lowest energy states, 
and a connecting exponential mesh of the same length con- 
taining the energetically higher states in the following way: 

_ / j ■ < j < N e 

kj -\N e + [ - l] J : N e < j < N t (B75) 

Here the symbol |_*J indicates the largest integer less than 
or equal to x, and £ <k 1 and 6 are adjustable parame- 
ters. The parameter 5 can be interpreted as the initial step 
size at the transition from the linear to the exponential mesh, 
i.e. fcjv e +i - kfj c ~ 5, which we have chosen as <5 = 2. The 
parameter £ is fixed by the requirement k^-\ = k max , leading 
to a transcendental equation for ^ which can be solved numer- 
ically. In our applications the complete mesh then consisted 
of N t = 80 points. We have generated an equivalent mesh for 
the sampling of the vibrational quantum number Kj, resulting 
in an 80 x 80 reduced kernel matrix 'K^ (Ej, B ), which has been 
used for the numerical solution of Eq. jB29t . Figure @ shows 
the reduced kernel matrix for a Vho = 300 kHz atom trap at a 
magnetic field strength of B = 158.1 G. The solutions f K of 
Eq. JB29> for the three lowest energetic three-body states ob- 
tained with the reduced kernel matrix of Fig.Qare illustrated 
inFig.E] 

6. Accuracy of the separable potential approach with respect 
to three-body energy spectra 

The range of validity of our exact solutions to the three- 
body Faddeev equations is limited by the accuracy of the inter- 
atomic potentials. Figure ^ clearly demonstrates that the si- 
multaneous adjustment of the separable potential of Appendix 
lAl to the scattering length of Eq. Q and to the formula @ 
for the binding energy of an alkali van der Waals molecule 
l42ll accurately describes both the measurements of Ref. I4r3l 
and their ab initio theoretical predictions |44|. This accuracy 
of the approach with respect to the energies of 85 Rb 2 Fesh- 
bach molecules persists over a wide range of magnetic field 
strengths, extending from the position of the zero energy res- 
onance of about 155 G up to 161 G, which is far beyond the 
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range of validity of the universal formula 0. Our separable 
potential approach also recovers the excited state energy spec- 
tra for a pair of 85 Rb atoms in a 300 kHz trap which we have 
determined using a microscopic binary interaction V{r) (see 
Fig.EJ. 

To demonstrate the validity of our separable potential ap- 
proach also in its applications to three-body energy spectra, 
we shall compare its predictions to the ab initio binding en- 
ergies of 4 He3 provided in Ref. fioll . As the low static elec- 
tric dipole polarizability of helium and the correspondingly 
small van der Waals coefficient of C(, = 1.461 a.u. 16511 do 
not allow us to accurately recover the helium dimer bind- 
ing energy from Eq. 0, we have used E^/ks = -1.313 mK 
(&b = 1.3806505 x 10~ 23 J/K is the Boltzmann constant) in 
addition to the scattering length of a — 190.7 a.u. to adjust the 
parameters A and cr of the separable potential of AppendixlAl 
These particular values of a and E\, correspond to the Tang, 
Toennies and Yiu (TTY) helium dimer potential ll65ll reported 
in Ref. Holl and determine the range parameter of the sepa- 
rable potential to be cr = 8.02 a.u.. This approach predicts 
the 4 He3 ground state energy to be E^/k^ = -96.9 mK and 
for the excited Efimov state we obtain E^ ep /k% = -2.09 mK. 
These predictions only slightly overestimate the exact ener- 
gies of £[ TY //fc B = - 126.4 mK and £J ty /A:b = -2.277 mK 
lUoll with relative deviations of 23 % and 8 %, respectively. 

To study the uncertainties of our approach to determine the 
separable interaction, we have also performed a different ad- 
justment of V sep based on the requirement that its effective 
range r^j p = {Acrl yfn)[\ - \fncrl{2a)] exactly recovers the 
effective range of r e g = 13.85 a.u. of the TTY potential in ad- 
dition to the exact scattering length. This determines the range 



parameter to be cr — 6.31 a.u.. The associated potential V sep 
yields E^ p /k B = -1.285 mK for the helium dimer, while it 
predicts the trimer energies to be E\ /&b = - 142.6 mK and 
E 2 /ku = -2.36 mK with relative deviations from the ab ini- 
tio results fioll of 12% and 4%, respectively. Both adjust- 
ments of the separable potential thus lead to similar degrees 
of accuracy, which are comparable to the accuracy of the adi- 
abatic hyper-spherical approach of Ref. @]. 

The small discrepancies between the predictions obtained 
from the different separable potentials and the ab initio calcu- 
lations of Ref. fioll indicate a remnant sensitivity of the trimer 
binding energies to properties of the microscopic binary inter- 
actions beyond those accounted for by the second order effec- 
tive range expansion of the two-body scattering phase shift. 
Such corrections can, in principle, be exactly included in our 
calculations by taking advantage of the universal properties 
of low energy three-body spectra discussed, e.g., in Ref. l66ll . 
This presupposes, however, that one of the trimer binding en- 
ergies is known either from experiment or from ab initio cal- 
culations. For instance, adjusting V sep in such a way that it 
exactly recovers the trimer ground state energy of Ref. HI Oil in 
addition to the binary scattering length yields cr = 6.79 a.u., 
which recovers the exact energy of the excited Efimov state re- 
ported in Ref. Jioll to a relative accuracy of 0.5 %. The studies 
of Ref. l67ll suggest, however, that the relatively small devi- 
ations between the ab initio and the different approximations 
employed in the separable potential approach may be physi- 
cally insignificant due to the remaining uncertainties of even 
the comparatively well known helium dimer interaction po- 
tential. 
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